Charge order and phase segregation in overdoped bilayer manganites 
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There have been recent reports of charge ordering around x — 0.5 in the bilayer manganites. At 
x — 0.5, there appears to be a coexistence region of layered A-type antiferromagnetc and charge 
order. There are also reports of orbital order in this region without any Jahn- Teller effect. Based 
on physical grounds, this region is investigated from a model that incorporates the two e g orbitals 
at each Mn site and a near-neighbour Coulomb repulsion, ft is shown that there indeed is both 
charge and orbital order close to the half-doped region coincident with a layered magnetic structure. 
Although the orbital order is known to drive the magnetic order, the layered magnetic structure is 
also favoured in this system by the lack of coherent transport across the planes and the reduced 
dimensionality of the lattice. The anisotropic hopping across the e g orbitals and the underlying 
layered structure largely determine the orbital arrangements in this region, while the charge order 
is primarily due to the long range interactions. 

PACS numbers: 75.47.Gk, 75.30.Et 



I. INTRODUCTION 



Bilayer manganites such as L^^aSri^xM^OV , 
the n = 2 member of the Ruddelsden-Popper se- 
ries {R, A) ri+ iMn n 0^ n +i (where R and A are rare- 
earth and alkaline-earth ions respectively), showi a 
distinct lack of symmetry in the nature of their 
ground states across x = 0.5 just as their 3D coun- 
terparts. These systems have started drawing at- 
tention after colossal magnetoresistance (CMR) has 
been reported in them 2 . These layered systems not 
only show large magnetoresistance and a sequence 
of magnetic phases 3 ^, they are very rich in their 
charge, magnetic and orbital structures. Most of the 
studies in the bilayer system so far focussed on the 
Mn 3+ richer region of a; < 0.5. The Mn 4+ — rich 
region (x > 0.5) is now explored carefully both 
experimentalljii*^ and theoreticallji^ and a succes- 
sion of magnetic phases C— > G has been con- 
firmed. The situation very close to and around x = 0.5 
is somewhat poorly understood yet. It emerged re- 
cently that there is a coexistence between charge or- 
dered and layered A-type spin ordered stated there, 
giving way to the C-type (or its polytypoi) magnetic 
order at larger Mn 4+ concentration. Both A-type 
(ferromagnetic layers coupled antiferromagnetically) 
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and C-type (ferromagnetic chains coupled antiferro- 
magnetically) phases have been found to be orbitally 
ordered. 

There are strong indications from several exper- 
imental groups that at and around x — 0.5, a 
charge ordered state coexists with an A-type anti- 
ferromagnetic phase. Both neutron scatteringifiiAi 
and muon spin rotation corroborate this view. A 
discontinuity in the muon precession frequency at 
x = 0.52 has been observed by Coldea et alJL Re- 
cently Wilkins et al.ji^ have reported that although 
there is clear evidence of both charge and orbital or- 
der in La2-2a;Sri + 2 ;c Mn207 at x = 0.5, there is no 
detectable signature for any Jahn- Teller (JT) distor- 
tion in their resonant X-ray diffraction study. JT dis- 
tortion of the MnC>6 octahedra begin to appear be- 
yond x — 0.55^. There is, though, a reasonably large 
(tetragonal) static distortion - the crystallographic c/a 
ratio reaches a minimumi around the same composi- 
tion, enhancing the AF super-exchange along c direc- 
tion. This is expected to favour the A-type spin order 
in that region. Unlike in the 3D manganites, there is 
no significant buckling of the bonds during distortion. 

It is now generally believed that the various mag- 
netic structures owe their origin to a large degree to 
the underlying orbital orderiii Models have been 
proposedMii^ for the manganites that incorporate the 
e g orbitals and the anisotropic hopping between them. 
The use of such models to the bilayer manganites (like 
La2-2a;Sri + 2 a ;Mn207 ) has only had limited success^ 
though, the A-phase is overestimated owing to the low 
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dimensionality of the lattice. A much improved mag- 
netic and orbital phase diagram was later obtained^ 
from a model that included the JT distortion as well 
as the proper choice of orbital degrees of freedom. The 
A-type AFM instability is indeed quite strong in the 
layered systemi, extending from x = 0.42 to 0.66 and 
at low temperatures. The FM region of the phase di- 
agram of some of the 3D manganites around x = 0.5 
is absent here, taken over by the A-type AFM phase. 
Having only two layers separated by large distance 
from the adjoining bilayer in the c-direction impedes 
a long range charge transport in that direction leading 
to stronger AF correlations along c-direction thereby 
contributing to an enhanced A-type instability. Be- 
yond x — 0.66, it is replaced by the C-type (and its 
poly type) order. The A- and C-phases are orbitally 
ordered and there is a connection between the pre- 
ferred orbital order and the observed magnetic order. 

However, the presence of charge order and its re- 
lation to the underlying magnetic phases have not 
been investigated in any of the previous theoretical 
studies^ in bilayer systems. The question of charge 
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As usual the charge and spin dynamics of the con- 
ventional DE modeli^ operate here too, with addi- 
tional degrees of freedom coming from the degen- 
erate e g orbitals (a, (3 stand for the e g orbitals) and 
hopping across them, which is determined by their 
symmetry. In the reasonably large Hund's coupling 
limit, which is relevant for the studies on mangan- 
ites, the DE mechanism implies that an electron can 
hop onto a site if core spins at that site (as well as 
the spin of any e g electron there) are parallel to its 
spin. Mobility reduces drastically with increasing J# 
if they are antiparallel. In Eqn.(l) Si and Si are the 
t2 g and e g spins at site i and Jh and Jaf are the 
Hund and super-exchange (SE) coupling respectively. 
Hint = U' J2iaa' r ^i<x"-i2cr' + V J2i ^3 contains the 
on-site inter-orbital and the near-neighbour Coulomb 
repulsion term. The intra-orbital term can be ignored 
as Jh in manganites is fairly large, preventing dou- 
ble occupancy in any given Mn 3d orbital. The ex- 
change interaction between two bilayers is known to 
be at least a 100 times weaker— than the intra-bilayer 
exchange. Two bilayers are also well-separated in an 



order and possible phase separation in the overdoped 
bilayer systems, therefore, still remains an open issue 
with increasing evidence in favour of such coexistence 
coming from the experiments. We attempt to address 
it in the following, starting from a model that incor- 
porates the essential physical attributes of this region. 



II. THE MODEL 

A quite general model for the bilayer manganites 
has been used£ to delineate the different orbital and 
magnetic structures for x > 0.5. It incorporates the 
degenerate e g manifold and the physics of double ex- 
change (DE) along with electron-electron interactions. 
This model can be adapted to investigate the charge 
order by including a longer range Coulomb term in the 
interaction part. Such a term has been knownSJiSi to 
give rise to coexisting charge ordered state in the 3D 
manganites at x = 0.5. 
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unit cell and intervened by the rare-earth ions. One 
can, therefore, use only one bilayer for a reasonable 
description of the system^. In general, in manganites, 
there is a strong JT coupling and that is included in 
the Hamiltonian. As we mentioned above, in the bi- 
layer systems (e.g. in La2-2xSri + 2iMn207 ), the JT 
coupling is weak and does not play a major rol o 6 i 12 
in the region x < 0.6. In the absence of electron- 
lattice coupling, the kinetic energy (KE) of electrons 
in the e g band and the Hund's coupling between ti g 
and e g spins compete with the antiferromagnctic SE 
interaction leading to a variety of magnetic and orbital 
structures. Considering that there are two e g orbitals, 
the nominal band filling is 

Typical values of the interaction and band param- 
eters for the bilayer systems are in the same range as 
in the 3D manganites. As in their 3D counterparts, 
the Hund coupling and Coulomb correlations are the 
largest scale of energjii& in the problem. We neglect 
the inter-orbital Coulomb term in the following dis- 
cussion of charge order. The inter-orbital Coulomb 
term does not have strong effect in the doping range 
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x less than or around 0.5 primarily because of the low 
filling ((1 — x)/4) making two orbitals at the same 
site less likely to be occupied simultaneously. As we 
see below, the magnetic order also prefers preferen- 
tial occupancies in one of the two e g orbitals in both 
A- and C-phases and the U' term only enhances that 
(discussed later). 

Treating the t2 g spins semiclassically 18 , the SE con- 
tribution to the ground state energy becomes E$e = 
j s 2 

A 2 (2cos8 xy + cos9 z ) where 9 xy and 9 Z are the angle 
between the near-neighbour t 2g spins in the xy plane 
and z direction respectively. In the Ferromagnetic 
state, xy — 9 Z = while in the A-phase 8 xy = 
and 9 Z — n. 

In the limit of infinite J#, the e g electron quan- 
tization axes at each site i is rotated in the local 
coordinate frame to make it parallel to Sj. This 
is accomplished quite easily by the spin- 1/2 rota- 
tion matrix exp(i^a z ) exp{i^fa y ) exp(—i^-a z ) oper- 
ating on a two component spinor. Allowing the core 
spins at each site to cant in the xz— plane and ne- 
glecting the phase term (the Berry phase) appear- 
ing from the transformation, it is straightfoward to 
showi£ that the effective hopping matrix elements are 
t xy — tcos{9 xy /2) and t z = tcos(9 z /2). In this level 
of approximation, the diagonalisation of the KE part 
of H reduces to solving the 2x2 matrix equation 
\\e a p — w6 a p\\ = for a system of spinless fermions. 
The matrix elements are obtained from the standard 
table of the overlap integrals^ involving d x 2_ y 2 and 
d 32 2_ r 2 orbitals on neighbouring sites of a square lat- 
tice and the dispersions are 

en = —2t xy (cosk x + cosky) 



£12 = *2i = j=2t xy (cosk x — cosky) 



£22 = -^t xy (cosk x + cosk y ) ~ -t z cosk z . 

Using the form of t xy and t z in terms of 9 xy and 9 Z in 
the Jh - > oo approximation, the two energy bands are 
obtained by the diagonalization of the above matrix. 
In the uncanted A- and C-phases, the dispersions for 
A- and C-phases become two- and one-dimensional 
respectively. In the pure FM phase, the DOS is three 
dimensional. However, in the bilayer system, owing 
to the absence of dispersion in the z direction, the 
DOS in both FM and A phases show two dimensional 
character as shown in Fig. 1. There being only two 
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FIG. 1: Density of states in pure FM and A phase. 



A;— points in the z— direction, the DOS for the C-phase 
has two delta functions centred at In all the 

calculations that follow, all energies are measured in 
terms of the overlap between the d 3z 2_ r 2 orbitals along 
z— direction t z = t, which has a typical value about 
0.25eV in manganites. 

III. CHARGE AND ORBITAL ORDERING 

In order to look for charge ordering in the bilayer 
manganites, we treat the nearest neighbour Coulomb 
interaction Vj^ < y > hihj in the mean field approxi- 
mation, < rii >= n + Coexp(iQ Ti) where Co is the 
charge order parameter and n is the average number 
of electrons per site. We take the usual staggered or- 
dering Q = (jr, 7r, 7r). The charge order parameter Co 
is then calculated self-consistently. The mean-field ap- 
proximation is known to work quite wellS^^ for the 
ground state properties in the manganites. In the infi- 
nite Jh limit, the electronic part of the Hamiltonian is 
a 4 x 4 matrix £ k Q/3 e ka /3cl a c kl3 - A £ k ct cl a £ k+Qa 
(where c represents the locally rotated e g electron op- 
erators described earlier) and this is diagonalized at 
each of the k points on a momentum grid. 

The bands obtained thereof are filled upto a chem- 
ical potential and the order parameter A = zVCq (z 
is the number of nearest neighbour) calculated along 
with the filling. The process is repeated until self- 
consistency is achieved as is customary in the mean- 
field theory. Charge order is indeed observed in a re- 
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FIG. 2: Charge order parameter as a function of V for 
three different angles 8 Z in the limit Jh —> oo. In the inset 
is shown the variation of charge order parameter with 9 Z 
for a fixed V = 0.7. 



gion x > 0.5 when V reaches a critical value, similar 
to the 3D manganites^Si. But unlike in the 3D case, 
we do not observe any F-phase coexisting with the 
CO in this region. As reported in previous work 8 ' 22 
A-phase instability is quite strong in the layered man- 
ganites owing to the 2D structure of its DOS and even 
for V — 0, there is no F-phase for x > 0.5i The 
A-CO coexistence region extends in a region above 
x = 0.5 (discussed later). In this limit of infinite Jh , 
the value of V for which the CO phase appears de- 
pends strongly on the canting angle 9 Z away from the 
A-phase (towards an F-phase). In Fig. 2 is shown the 
the dependence of the CO order parameter on V at 
three different angles (6 Z = 180 deg being the pure A- 
phase). The inset shows how Co varies with canting 
at a typical V = 0.7. In the infinite Hund's coupling 
limit, the only way an electron gains KE in and AF 
configuration of background spins is via the canting 
of them 2 ^, thereby generating an effective non-zero 
hopping across. This has been observed^ in the 3D 
manganites as well close to its G-C phase boundary 
(for x ~ 1). As we observe here, a CO phase does not 
coexist with an F phase (which is energetically un- 
favourable to an A-phase in this region) and therefore 
the increased canting requires larger V to bring about 
the charge ordering. At a fixed V, therefore, the CO 
order parameter reduces with increased canting away 
from A-phase. 



The appropriate limit for any of the manganites 
is, of course, a large but finite value of the Hund's 
coupling as is generally believed. In this limit, we 
can treat the core spins semi-classically again and 
the spin degrees of freedom for the e g electrons are 
reintroduced. In order to take care of this, the sec- 
ond term in the Hamiltonian (1) would have to be 
treated now in various possible ground states. We 
choose Si = S exp(iq.ri) to represent a homogeneous 
spin configuration, where q determines different spin 
arrangements for the tig spins 20 . The second term 

in eqn.(l) becomes — JhSq X/kaer < - rc kQ<T C i <; +qacr- The 
near-neighbour Coulomb term is treated in the mean- 
field as above. In this semi-classical approximation for 
the ti g spins the Hamiltonian (1) reduces, then, to a 
16 x 16 matrbi^ii at each k point. 

This mean-field Hamiltonian is again diagonalised 
at each k-point on the momentum grid. The ground 
state energy is calculated for different magnetic struc- 
tures. The CO order parameter is also determined 
self-consistently. Four different magnetic structures 
are relevant for the experimental phase diagram (with 
q values in the parentheses), A-type (0,0, w), C-type 
(-7T, 7r, 0) - the usual C-phase with FM chains along c- 
direction, and the 3D AFM G-type (7T, 7t, 7t). There 
is also the C'-type polytype structure^, with q = 
(-7T, 0, tt) This is same as a C-type, only that its FM 
ordering is along y-direction as reported by Ling et al. 
The magnetic structure with minimum ground state 
energy and the CO order parameter are determined 
for each set of parameters {x, Jh, Jaf) for the range 
of doping (0.5 < x < 1). 

Charge order is observed in the region 0.67 > x > 
0.5. Fig. 3 shows the CO order parameter as a func- 
tion of V for x — 0.5. The underlying spin order in 
this region is found to be A-type and there is a A-CO 
to A transition at around x = 0.67 for V = 0.7. There 
is a jump in the CO order parameter as a function of 
x, shown in Fig. 4, indicating a first order transition. 
The coexistence region of CO and A is a manifestation 
of this first order transition. Based on the observed 
variation of these different charge and spin ordered 
states, we obtained a phase diagram in the x — Jh 
plane for V = 0.7. 

The observed phase diagram, shown in Fig. 5, is 
markedly similar to the one obtained by Ling and 
coworkers. In addition it also has the charge ordered 
phase coexisting with the A-phase in the region close 
to x = 0.5. As reported by Coldea et al.^, there is 
indeed a region of A-CO coexistence between x = 0.5 
and x = 0.65. We note that we have also tried a 
charge order with wave vector Q = (77, 7r, 0) instead of 
the fully staggered one and observed that energetically 
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FIG. 3: Charge order parameter as a function of V at 
x = 0.5. 



FIG. 5: Phase diagram, in the JhSo versus filling x plane. 
The terms A, C', G, CO etc. are explained in the text. 
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FIG. 4: Variation of charge order as x is changed from 
0.5 for a fixed V = 0.7. The jump in the order parameter 
indicates the A-CO to A boundary 



it is very close to the staggered one. This is primarily 
due to the absence of coherent charge transport in the 
c-direction making a charge uniform state in that di- 
rection nearly degenerate with a charge ordered one. 
We note in passing that we did not observe the well 
known CE phased* in our calculation even at x = 0.5. 
This is similar to the findings in the V = case stud- 



ied earlienZA The region of no long range order just 
above the A-CO coexistence reported by both Ling 
et al. and Coldea et al. is beyond the scope of the 
treatments here. It is at this point useful to note that 
a region of two phase coexistence is not very stable 
against long range Coulomb interactions and in real 
systems one would probably observe a microscopically 
phase separated mixture of one phase in another. The 
region of no LRO could well be such a region, around 
the A-CO to A transition and with the competing C 
phase energetically very close by. 

In the absence of CO the phases A and C are both 
orbitally ordered^ and we look for this in the pres- 
ence of the near-neighbour Coulomb term also. We 
do indeed find very similar orbital occupancies here 
too. There is a predominant occupancy in the d^ z -2_ r i 
orbital over the d x 2_ y 2 orbital indicating orbital or- 
dering. The occupancies in the two orbitals must add 
up to the total electron density (-^-j 2 )- The effect of 
charge ordering is seen to be small, shown in Fig. 6, 
reducing the orbital order (measured by the imbal- 
ance in the orbital occupancies) in the A-CO phase 
from its value in the absence of CO. There is no ef- 
fect in orbital occupancies at the A-CO to A bound- 
ary, they smoothly continue into the pure A-phase. 
The presence of orbital order here is not contingent 
upon the underlying lattice distortion or JT effects. 
This conforms to the recent observation of Wilkins et 
alH that even without any JT distortion, there is a 
pronounced orbital order. In the region of C'-phase, 
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FIG. 6: Orbital occupancies in the A-CO and A-phases. 
The upper curves are for d 3z 2_ r 2 while the lower ones are 
for the d x 2_ y 2 orbitals. The dotted lines are for compari- 
son with V = situation where no CO is present. 



the orbital order is identical (predominantly d 3y 2_ r 2 
over d z 2_ x 2) to the one reported earlierifor the undis- 
torted bilayer system and are not shown here. We 
have looked into the effect of an inter-orbital repul- 



sion term on the phase diagram and charge ordering 
and except for a slight enhancement of the orbital or- 
der, found it to have little effect on the phase diagram. 
The intra-orbital Coulomb term is known to have very 
little effectAS, in the mean-field theory in the region of 
large J# and we haven't considered it in the above. 



IV. CONCLUSION 

Motivated by the observation of charge ordering in 
the region x > 0.5 in the bilayer systems, we have 
investigated it from a model that has been quite use- 
ful in understanding the overdoped manganites. The 
charge order has been found in our calculation in the 
same region where it has been seen experimentally. 
There is also a region of phase separation abutting 
the experimentally observed no LRO region. We ob- 
serve orbital order even in the absence of JT distortion 
as reported recently in the same region of doping. It 
would be interesting to see more experiments on the 
region of no LRO in the phase diagram. The possi- 
bility of canting of spins is yet not ruled out and we 
have seen trends for it for very large Hund's coupling 
in our calculations. 
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